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Abstract 

Ground-state properties of fermionic mixtures confined in a one-dimensional 
optical lattice are studied numerically within the spinless Falicov-Kimball 
model with a harmonic trap. A number of remarkable results are found, (i) 
At low particle filling the system exhibits the phase separation with heavy 
atoms in the center of the trap and light atoms in the surrounding regions, 
(ii) Mott-insulating phases always coexist with metallic phases, (iii) Atomic- 
density waves are observed in the insulating regions for all particle fillings near 
half-filled lattice case, (iv) The variance of the local density exhibits the uni- 
versal behavior (independent of the particle filling, the Coulomb interaction 
and the strength of a confining potential) over the whole region of the local 
density values. 

PACS numbers: 05.30.Fk, 71.10.-w, 71.30.+h 
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1 Introduction 



The ability to confine ultracold Bose and Fermi gases inside artificial crystals gener- 
ated by standing- wave laser light fields, i.e., optical lattices, offers the possibility to 
create ideally clean and highly tunable strongly interacting quantum many-body sys- 
tems pp . The low-energy properties of these systems can be described using models 
borrowed from condensed matter systems [2] whose parameters and dimensionality 
can be controlled with high precision. This fact opens new routes for understand- 
ing the physics of long-standing problems of strongly correlated systems, like the 
phase separation, metal-insulator transitions, superconductivity, etc. On the other 
hand, the quadratic confining potential, present in addition to the regular "lattice" 
potential leads to a number of fundamentally new phenomena. For example, it was 
shown [3, 4j that in the presence of a confining potential the Mott-insulating phase is 
restricted to the domain that coexists with a metallic phase, in contrast to the global 
character typical of solid-state systems. Moreover, mixtures of two-component atoms 
with different masses (e.g., 6 Li and 87 Kb) introduce an additional parameter, namely, 
the difference between the hopping amplitudes associated with each species of atoms 
in the optical lattice. This may affect the stability of the possible quantum phases 
or even induce new ones [5J El U\- For these reasons the mixtures of ultracold atoms 
in optical lattices belong to the most intensively studied subjects of contemporary 
experimental and theoretical physics. 

In this paper we investigate the ground-state properties of mixtures of two species 
of fermionic atoms with strongly different masses in a harmonic potential. Since the 
tunneling rate decreases exponentially with the square root of mass of the atom, 
and particles of different species on the optical lattice interact only through the on- 
site interaction, this leads naturally to the Falicov-Kimball Hamiltonian [S] with a 
confining potential [3] 

H = +Uj24d i ftfi + (£)V£(*-§) (dfd, + f? fr) (1) 

where ftifi) and df(dj) are the creation (annihilation) operators of heavy (/) and 
light (d) particles at lattice site i. The number of lattice sites is L and is selected 
so that all the fermions are confined in the trap. We denote the total number of 
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fermions in the trap as iV and consider equal number of heavy (JV/ = ftfi) an d 
light (Nd = J2i dfdi) atoms. 

The first term of Eq. 1 is the kinetic energy corresponding to quantum- mechanical 
hopping of the light d atoms between the nearest-neighbor sites i and j. These 
intersite hopping transitions are described by the matrix elements ty, which are 
— td if i and j are the nearest neighbors and zero otherwise (in the following all 
parameters are measured in units of td). The second term represents the on-site 
Coulomb interaction (U > 0) between the light and heavy atoms. The last term is the 
energy of light and heavy atoms in the harmonic trapping potential. In accordance 
with similar studies for the asymmetric Hubbard model [9J, we consider here the 
same trapping potential for both species of atoms. 

Since in this spinless version of the Falicov-Kimball model with a confining po- 
tential the /-heavy atom occupation number ff fi of each site i commutes with the 
Hamiltonian (1), the /-heavy atom occupation number is a good quantum number, 
taking only two values: Wi = 1 or 0, according to whether or not the site i is occupied 
by the heavy atom. Therefore, the Hamiltonian (1) can be written as 

h = £ Mto + (£)V£(;- ^fw h (2) 

where h ij (w) = ^ + {uwi + (~) V (i - §) ^ <%. 

Thus for a given configuration of heavy atoms w = {wi,w 2 ■ ■ ■ wi} defined on the 
one-dimensional lattice, the Hamiltonian (2) is the second-quantized version of the 
single-particle Hamiltonian h(w), so the investigation of the model (2) is reduced to 
the investigation of the spectrum of h for different configurations of heavy atoms. 

It is well known that in the absence of harmonic confinement (V = 0) the 
one-dimensional Falicov-Kimball model exhibits a rich spectrum of solutions includ- 
ing phase separated (U small) as well as most homogeneous distributions of heavy 
atoms [10]. However, due to the confining potential the lattice sites become, inequiv- 
alent, and thus, it is of fundamental importance to analyse the interplay between the 
on-site Coulomb interaction and the confining potential. 

To describe the system at nonzero V we have calculated various local quantities, 
like the local density of heavy atoms (n{ = (/j + /i)), the local density of light atoms 
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{nf = (dfdi)), the total site occupation (n* = {ft fi + dfdi)), the variance of the local 
density (Aj = {{ft ft + dfdi) 2 ) — {ft fi + dfdi) 2 ) and the local double occupation 
{Di = {dfdiftfi)), as functions of total number of the confined atoms, the Coulomb 
interaction U and the confining potential V. The ground states are calculated by a 
well-controlled numerical method that we have elaborated [TT] for a description of the 
conventional Falicov-Kimall model (V = 0). Later, the method was successfully used 
for various generalizations of the Falicov-Kimball model and different physical prob- 
lems [12]. Its generalization on systems with a harmonic potential is straightforward 
too. 

2 Results and discussion 

In Fig. 1 we present results of our numerical calculations for n{ ,nf,rii, Di and Aj 
obtained on the one- dimensional cluster of L = 120 sites at U = 4, V = 4 and 
different fillings. We have added also the profiles of the local compresibility that 
has been proposed by Rigol et al. |3j as a local order parameter to characterize the 
Mott-insulator regions. This quantity is defined as [3] 

K l = Z! Xi,i+j , (3) 

\j\<l(U) 

where 

Xi,j = {niUj) - {ni){nf) (4) 

is the density-density correlation function and l(U) ~ b£(U), with £(U) the correla- 
tion length of Xi,j i n the unconfined system at half-filling for a given U. The factor 
b is chosen within a range where k 1 becomes qualitatively insensitive to its precise 
value [4J. The insulating regions are then characterized by k\ = 0. For the values of 
U used here we usually have b ~ 4-8 with £(£/) ~ 1. 

The most interesting result obtained at low particle fillings (we note that Nf = 
Nd = N/2) is the observation of the complete phase separation with heavy atoms 
in the center of the trap (n{ = l,nf ~ 0) and light atoms in the surrounding 
regions (n{ = 0,nf < 1). Such a behaviour is found for all particle fillings from 
Nf = 1 to some critical value NS that rapidly decreases with increasing V and 
is almost independent of the local Coulomb interaction U. In the regions where 
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n{ = (0 < nf < 1) the variance of the local density Aj and the local compresibility 
k\ are finite (the metallic phase), while in the middle of the trap where n{ = 1 (nf 
exponentially decreases in this region) both Aj and k\ are equal to zero (the insulating 
phase). Thus in accordance with results obtained for the Hubbard model [3J (the 
hopping probabilities are same for both types of atoms) we have found that also 
in the Falicov-Kimball model insulating domains coexist with metallic regions, such 
that global quantities are not appropriate to describe the system. 

At higher particle fillings the situation is more complex. Above the critical filling 
Nj (Nf = 24, for U = 4, V = 4) the connected cluster by heavy atoms occupied 
sites (in the center of the trap) splits on smaller clusters, usually of the same size, 
separated by the empty site. As Nf increases the size of clusters decreases from Nf 
to 1. Of course, the redistribution of heavy atoms has dramatic consequences on the 
distribution of light atoms. Now, the light atoms occupy preferably the empty sites 
in the middle of the trap (n, ~ 1) what leads to the atomic-density waves in nf and 
ni profiles. In the region where ~ 1 the variance of the local density is finite but 
smaller than in surrounding metallic regions indicating [3j H] the presence of Mott- 
insulating phase in the center of the trap. This conjecture supports the behaviour of 
the local compresibility k\ that is equal to zero in the corresponding region. 

Increasing the number of particles up to Nf = 52, the connected cluster of heavy 
atoms (n{ = 1) starts to form in the middle of the trap. In this region both the 
variance of the local density and the local compresibility are finite what indicates the 
presence of metallic phase in the center of the trap. Upon adding more fermions, this 
new metallic phase widens spatially, while the Mott-insulating regions of the atomic 
density waves are pushed to the borders and completely disappear at Nf = 90. 
At this filling practically the whole region of the trap, except the edges (where 
n{ — 0, nf ~ 1), is metallic. With still a higher filling this metallic phase is further 
stabilized, but at some critical filling (Nf = 95) a new insulating phase ("a band 
insulator") starts to develop in the center of the trap (n{ = l,nf = 1). This trend 
holds also for the highest particle fillings, the width of the band-insulating phase 
increases and the surrounding metallic regions are gradually suppressed. 

To reveal the role of the Coulomb interaction U and the confining potential V on 
a formation of metallic and insulating domains, similar calculations have been also 
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performed for various combinations of U and V. The results of numerical calculations 
for the site occupation rt, as a function of i are displayed in Fig. 2. One can see that 
the Coulomb interaction U and the confining potential V exhibit precisely opposite 
effects on the stability of metallic and insulating domains. Indeed, with increasing V 
(at fixed U) the insulating domain (n$ ~ 1) is suppressed and the metallic domains are 
stabilized, while with increasing U (at fixed V) the metallic domains are suppressed 
and the insulating domain is stabilized. 

Since the Falicov-Kimball model can be considered as a simplification of the 
Hubbard model (only one kind of particles, say with spin up can hop) it is interesting 
to compare results obtained in these two different limits. Such a comparison (see 
Ref. 3 and Ref. 4) reveals obvious differences in behavior of these models in the 
confining potential. For example, the ground-state of the Hubbard model (for V ^ 0) 
is always metallic at low particle fillings, while the metallic regions coexist with the 
insulating region in the ground state of the Falicov-Kimball model. Moreover, the 
local density profiles exhibit obvious oscillations for the Falicov-Kimball model, while 
no sign of such oscillations has been observed for the repulsive Hubbard model (with 
the exception of Hartree-Fock [4 J and variational studies [13] , that were not confirmed 
by a projector Monte Carlo simulations, however [4]). 

To exclude the possibility that oscillations are a consequence of a finite size of 
clusters used in our numerical calculations, we have performed an exhaustive finite- 
size scaling analysis on finite clusters up to L = 480 sites for all particle fillings 
from Fig. 1. This analysis showed that the ground states found for L = 120 hold 
also on clusters of L — 240 and L = 480 sites and thus they can be satisfactorily 
extrapolated on much larger clusters. In Fig. 3 we present numerical results for 
extrapolated ground states obtained for a cluster of L = 6000 sites and the same 
values of rif = Nf/L, V and U as used in Fig. 1. These results clearly demonstrate 
that increasing L suppresses the atomic-density oscillations in the metallic phase, 
but stabilizes the atomic-density waves in the insulating phase. 

Since in the real experiments with ultracold atoms the hopping matrix elements 
(of heavy atoms) between the nearest-neighbor sites (tf) are not strictly equal to 
zero, it is necessary to examine the stability of our solutions obtained for tf = 
against the finite values of tf. For this reason we have performed exhaustive exact- 
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diagonalization studies of the asymmetric Hubbard model [H] (tf 7^ 0) in a confining 
potential for a wide range of model parameters (Nf,tf,U and V) on finite clusters 
up to L = 12. The representative examples of local density profiles are displayed in 
Fig. 4 for several different values of tf. They clearly show that results obtained for 
tf = remain stable also at finite (small) tf. Thus the Falicov-Kimball model with 
a confining potential can be used satisfactorily to model the ground-state behavior 
of mixtures of two-component fermionic atoms (with strongly different masses) in a 
harmonic potential. 

Finally, we have also calculated the variance of the local density as a function 
of the local density rii. It is known from the study of the one-dimensional Hubbard 
model in a harmonic potential |[3 J that this quantity shows universality with respect 
to the confining potential for rii — > 1. It should be noted that the universal behavior 
was observed only for strong Coulomb interactions U, where systems have a Mott- 
insulating phase at n = 1. In Fig. 5 we present numerical results for the variance 
of the local density obtained within the Falicov-Kimball model with a harmonic 
potential as a function of the local density calculated for various particle fillings and 
various values of U and V. One can see that all numerical data for the variance of 
the local density collapse on the same curve (given by A = —n(n — 1) for n < 1 
and A = — (n — l)(n — 2) for n > 1). Thus in contrast to the similar studies on the 
Hubbard model [3] we have found that the variance of the local density of systems 
described by Falicov-Kimball model with a harmonic potential exhibits universality 
not only with respect to the confining potential but also with respect to the Coulomb 
interaction U. In addition, we have found that the universality hold not only for 
rii — > 1, but for all values of n^. A similar universal behaviour has been observed 
also for other local quantities. For example, the local double occupation Di = for 
rii < 1 an d Di = rii — 1 for n.j > 1, independent of values of U, V, Nf and L. Also the 
universal behavior is found for the local compresibility k 1 when rij — > 1. However, 
the value of the critical exponent for the Falicov-Kimball model is equal to 1, unlike 
the nontrivial value 0.68-0.78 found for the Hubbard model [3J. 

In summary, we have studied the ground-state properties of fermionic mixtures 
with mass imbalance in a one- dimensional optical lattice within the spinless Falicov- 
Kimball model with a harmonic potential. We have found that the system exhibits 
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the phase separation at low particle fillings. In this case the heavy atoms occupy the 
center of the trap while the light atoms are localized in the surrounding (metallic) 
regions. At higher fillings we have observed a formation of Mott-insulating domains 
with atomic-density waves. In all cases Mott-insulating phases coexist with metal- 
lic phases. One of the most interesting results is, however, the observation of the 
universal behavior of the variance of the local density (independent of the particle 
filling, the Coulomb interaction and the strength of the confining potential) over the 
whole region of the local density values. 

This work was supported by Slovak Grant Agency VEGA under Grant No. 2/7057/27 
and Slovak Research and Development Agency (APVV) under Grant LPP-0047-06. 
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Figure Caption 



Fig. 1. Profiles of various local quantities along the trap calculated for V = 
4, U = 4, L = 120 and four different values of particle fillings. Different panels (from 
top to bottom) correspond to: the local density of heavy atoms n{, the local density 
of light atoms nf, the total site occupation n» = n{ + nf, the local double occupation 
Di, the variance of the local density Aj and the local compresibility k\. 

Fig. 2. The site occupation function of site position % calculated for 

different V (the first column) and different U (the second column). 

Fig. 3. The site occupation function of site position % calculated for 

extrapolated ground states for different particle fillings. The parameters involved 
are V — 4, U = 4 and L = 6000. The insets show details of rii profiles for the region 
of atomic-density waves. 

Fig. 4. The exact-diagonalization results (the asymmetric Hubbard model with 
a harmonic trap) for the site occupation as a function of site position % calculated 
for different values of tf, Nf, U and V on a finite cluster of L = 12 sites. 

Fig 5. The variance of the local density A as a function of the local density 
n calculated for different values of U, V and Nf. The numerical data are fitted by 
A = — n(n — 1) for n < 1 and A = — (n — l)(n — 2) for n > 1. 
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